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Abstract 

Early growth of density fluctuations of nuclear matter in spinodal region is investigated employing 
the stochastic mean-field approach. In contrast to the earlier treatments in which only collective 
modes were included in the calculations, in the present work non-collective modes are also included, 
thus providing a complete treatment of the density correlation functions. Calculations are carried 
out for symmetric matter in non-relativistic framework using a semi-classical approximation. 
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I. INTRODUCTION 


Much work has been done to analyze multi-fragmentation processes observed in heavy-ion 
collisions. These processes exhibit many features compatible with liquid-gas phase transition 
of nuclear matter. Experimental measurements of charge correlations between produced 
fragments suggest that the induced phase transformation is due to the amplification of 
density fluctuations in the spinodal region. Microscopic description based on mean-field 
approximation can provide only a limited information about the spinodal dynamics. In 
the mean-field framework, it is only possible to determine the boundary of the dynamical 
instabilities and the growth rates of the most-unstable collective modes. Growth of density 
fluctuations and large amplitude dynamical evolution of the system require a microscopic 
framework beyond the mean-field that incorporates the dynamics of density fluctuations 

[1-5]. 

Considerable effort has been spent to improve the standard mean-field approximation 
by incorporating dynamics of density fluctuations [6-9]. In a recently proposed stochastic 
mean-held approach (SMF), the effect of quantal and thermal fluctuations in the initial 
state is incorporated into the description in a stochastic manner [10]. The standard mean- 
held approximation is a deterministic description in the sense that a well- defined initial 
state gives rise to a unique hnal state. On the other hand, in the SMF approach, starting 
from a proper distribution of initial states, specihed by quantal and thermal fluctuations, an 
ensemble of mean-held trajectories are evolved each with its own self-consistent Hamiltonian. 
A number of recent publications have provided strong validation of the SMF approach for 
describing dynamics of density fluctuations at low energies [11-16]. 

Employing the SMF approach in non-relativistic framework [17], as well as in relativistic 
framework based on Walecka-type effective held theory [18], we recently carried out a num¬ 
ber of investigations of the early development of spinodal instabilities and baryon density 
correlation functions in nuclear matter in the semi-classical and quantal descriptions [19-22], 
In the linear response regime, linearizing the equation of motion around a suitable initial 
state inside the spinodal zone and employing the method of one-sided Fourier transform, it is 
possible to carry out nearly analytical treatment of the baryon density correlation functions. 
In these investigations only unstable collective modes were included in the calculations. As 
pointed out in Ref. [23], for a complete treatment of initial growth of density fluctuations it 
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is necessary to include non-collective modes as well. In the present investigation, we under¬ 
take calculations of correlation functions of baryon density fluctuations including collective 
as well as non-collective modes. 

In Section 2, we present nearly analytical description of the density correlation func¬ 
tions, including collective and non-collective modes in the linear response framework, in the 
semi-classical limit. In Section 3, we present calculations at different initial densities and 
temperatures for symmetric matter. Conclusions are given in Section 4. 


II. DENSITY CORRELATION FUNCTIONS IN SPINODAL REGION 


The SMF approach goes beyond the standard mean-held theory by including quantal and 
thermal fluctuations in the initial state. An ensemble of single-particle density matrices are 
generated by incorporating the initial fluctuations. It is possible to formulate the approach in 
relativistic and non-rclativistic frameworks. Here we consider the non-relativistic approach. 
The single-particle density matrix of each event (indicated by the event label A), p x (r,r',t), 
is evolved through to its own mean-held according to 


ih 



ha(t),Pa(t) 


( 1 ) 


where the label a — p i,n t , n j. indicate the isospin-spin quantum numbers, and h x (t) 
is the mean-held Hamiltonian for the event. For the description of early growth of density 
fluctuations around an initial state p a , it is sufficient to linearize the equation of motion 
around this initial state. Subsequently small fluctuations, Sp x (t ) = p x (t) — p a , in the density 
matrix are determined by 


ih^Spait) = [ha,Spa(t)] + [6U*(t),p. 


( 2 ) 


where h a denotes the mean held Hamiltonian at the initial state and SU X is the huctuating 
part of the mean-held for the event. We are interested in investigating the development of 
spinodal instabilities in nuclear matter around a uniform initial state p a . In this case, by 
employing a plane wave representation, we can carry out an almost analytical treatment 
of the density correlation function. In the plane wave basis, the linear response equations 
become 
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(3) 


ih§- t <Pi\5p x a {t)\p-2 > 

= [£o(pi) - £a(P2)\ < Pl\Spa(t)\p2 > ~ [fa(Pl) ~ fa (&)] < Pl\$U X (t)\p 2 >, 
here f a (p) denotes Fermi-Dirac factors at finite temperature. In the current work, we con¬ 
sider symmetric nuclear matter and assume that the mean-field potential depends only on 
the local nucleon density. Consequently, small fluctuations in the mean-held can be writ¬ 
ten as SU x (r,t) = ( dU/dp)o5p x (r,t ), where 5p x (r,t ) denotes the small fluctuations in total 
nucleon density. Also, because of relative simplicity, we carry out calculations in the semi- 
classical limit. The Fourier transform of nucleon density fluctuations, Sp x (k,t) is related to 
the fluctuations of the density matrix according to 

Sp x (k, t) = J2 [ L < P + h ^/ 2 \ 6 Pa\P ~ kik/2 >, (4) 

a 2 ( 27 \Ti) 

where summation a is over the spin-isospin quantum numbers. The method of one-sided 
Fourier transform is very useful for solving Eq. [3] [24,25]. We introduce the one-sided Fourier 
transform of the local density in time as 

5p x (k,uj) = / dte lut 5p x (k : t). (5) 

Jo 

This yield an algebraic equation for the Fourier transform of the local density, which can be 
solved to obtain 


5p x (k,uj) = — -4 — -G x (k,uj), ( 6 ) 

£{k,UJ) 

where e(k,uj) = 1 + F 0 x(k,co) is the susceptibility, with F 0 = ( dU/dp)o as the zeroth order 
Landau parameter, and x(k,uj) denotes the Lindhard function. The semi-classical expression 
of the Lindhard function is given by 


The quantity G x (k,co) in 


iT \ A [ d3 P v-k df Q 

X{k,uj) = -4 / ^ r - -w-. 

J (27 Th) s y . p — cj as 

6 is determined by the initial conditions 


G\k,u) = Y. 

a 


d 3 p < p + hk/2\Sp x (0)\p — hk/2 > 
(2t rh) 3 v-k-u 


(7) 

( 8 ) 


which acts as a source for developing fluctuations. According to the basic postulate of the 
SMF approach, elements of the initial density matrix are uncorrelated Gaussian random 
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numbers with zero mean values and with well defined variances. In the semi-classical limit 
their variances are given by 


< p + hk/2\5p x (0)\p — Tik/2 >< p 1 — Tik'/2\5pl(0)\f? + Tik' /2 > 
= 5 ab (2nh) 3 5(p- p , )(27r) 3 5( y k - k') [f(p) (1 - f(p))\. 


(9) 


Here and below the bar over the expressions indicates the ensemble average. If we ignore the 
Coulomb potential, this has the same form for neutrons and protons in symmetric nuclear 


matter, therefore we omit the label a in the Fermi-Dirac factors. The factor 5 ab reflects 


the assumption that local density fluctuations in spin-isospin modes are uncorrelated in the 
initial state. 

According to the method of one-sided Fourier transform, time evolution of the density 
fluctuations is determined by the inverse transformation in time, which is expressed as a 
contour integral in the complex cu-plane as 



( 10 ) 


where integration path passes above all singularities of the integrand, as shown by line C\ 


in Fig.l. We can calculate this integral by employing the residue theorem and closing the 
contour in a suitable manner. For the calculation of this integral, we need to investigate 
analytical structure of the integrand. In the spinodal zone, there are collective poles de¬ 


termined by imaginary solutions of the dispersion relation e(k,uj) — 0 —> cu — Air*,. The 


collective poles play an important role in early growth of density fluctuations. However, 
as pointed out by Bozek [23], collective poles alone do not give the full description of the 
growth of instabilities. In fact, collective pole contributions alone do not even satisfy the 
initial conditions. By calculating the angular integral in the Lindhard function given by 
EqJTJ it is possible to see that, there is a cut singularity of the integrand in EqJlO] along 
real axis in the complex w-plane. In addition, G x (k,co) may also have singularities on the 
real axis. These non-collective poles of the integrand in EqjlO] can be seen more clearly in 
discrete description of the plane wave representation. If the system is placed in a finite box, 
plane waves are characterized by discrete values of momentum. In this case as illustrated 
in [26], dispersion relation has both imaginary collective poles and non-collective poles on 
the real axis. In the continuous limit, these non-collective poles appear as a cut singularity. 
Since the integrand in complex cu-plane is multivalued, the entire real cu-axis is a branch cut. 
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In order to calculate the integral in EqJTOj we choose the contour C , as shown in Fig. 1. 
We exclude the real cu-axis by drawing the contour from +oo to the origin just above the 
real ca-axis, and after jumping from the first Riemann surface to the second at the origin, 
we draw the contour just below the real cu-axis from origin to +oo . Contour is continued 
with a large semi-circle and it is completed by jumping from the second Riemann surface to 
the first one at origin,as shown in Fig.l. As a result, the integral can be expressed as 



FIG. 1. (Color online) The contour chosen to calculate the integral in Eq. 10 


the collective poles of the integrand, and the dotted line along the entire real 
cut of the integrand. 


Dots indicate 
w-axis is the branch 


5p X (k,t) = 3p X ole (k, t) + 5p X ut (k,t), 
where the Erst term is the pole contribution 




G x (k,±iT k ) 




de{k,u) / duj\ u =±ir k 


and the cut contribution is given by 


S Pcut( k ,*) = ~i 


r+°° du 
I- oo 27r 


G x (k, oj + ie) G x (k,u — ie ) 


—iu)t 


(a) 


( 12 ) 


(13) 


e(k,u + ie) e(k,u — ie) 

Local nucleon density fluctuations 5p x (r',t) are determined by the Fourier transform of 
5p x (k , t) as 


5p x (r,t) = 


d 3 k 
(2t r) 3 


e ik ^5p x (k,t). 


(14) 
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In order to investigate the development of the spinodal instabilities a very useful quantity 
is the equal time correlation function of local density fluctuations 

a(\r-f*\,t) = 5p x (r,t)5p x (r',t ) = j a(k,t). (15) 

The spectral intensity a(k, t) of the correlation function is determined in terms of the variance 
of the Fourier transform of density fluctuations according to 


a(k,t)(27r) 3 5(k — k') = 5p x (k,t)5p x (—k',t). (16) 

We can evaluate the spectral intensity a(k,t) by evaluating the ensemble average using the 
Eqs. (12) and (13) for pole §p x ole (k } t ) and cut 5p* ut (k,t) part of the Fourier transform of 
density fluctuations and the Eq. (J9]) for the initial fluctuations. As a result, the spectral 
intensity becomes 


a(k,t) = <J pp (k,t) + 2a pc (k,t) + a cc (k,t), (17) 

where the first and last terms are due to pole and cut parts of the spectral intensity and the 
middle term denotes the cross contribution. The pole part is 


(Tpp(/c, t) 


E 4 


de(k,uj)/duj 


(e +2Fkt + e - 


- uj=irk 


2 r fe t\ 

2 E_ 


) | 

d£(k,u)/du 

uj=iFk ^ 


(18) 


where 


E± = 4 


(1 “ m) 


(y-k) 2 ±n 

(v-ky + ri 


(19) 


The pole contributions are the same as what we found in our previous investigation [17]. 
The cut-cut contribution has four terms 


crcc(k,t) = A^(k,t) + A^(k,t) + A±(k,t) + A±(k,t), 


( 20 ) 


with 


A=f(k, t ) — 


r+°° du r+°° du' cp(u =f irj) + <j)(u' =F irj) 

’ —oo 2vr J —OO 2vr £ (k, u =F iri)e{%, u' =F irj) u + u> ^ 2 il i' 




and 


^ (k t) = - / + °° — [ + °° duj ' ^ T ir} ^ + y ± ir] ^ 1 c r i{u,+u,)t 

T ’ 7-oc 2 tt J -00 2tt £(fc ; (j =p irj)£^k, u' ± ip) w + W 


( 21 ) 


( 22 ) 
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In these expressions p is an infinitesimal positive number, and quantity <p{u j + ir]) is defined 
by 3 

<l>{v + iv) = 4 f f(p) (1 - /(/)) 1 -(23) 

j (2vr ny y.k-tu-n-j 

In the double integral of A^(k,t), there are principle value and delta function contributions 
which are identified using the identity 1/ (cu + to' + irj) = P( 1/uj + a/) ± ip5(cu + a/). The 
integrand of A T (k, t ), in contrast to its appearance, is a well behaved function, because when 
u / = — u), the nominator is also zero therefore the ratio (0(a; + ir /) + 0(a/ ± ir /)) / (u> + uo') 
becomes finite. 

Mixed term in the spectral intensity of Eq. (JTtJ) has the following form 
cr pc {k, t) = ( k, t ) + B t ( k, t) + B±(k,t) + B±(k, t ), 

with 

B (% t \ = ±i e ^ kt [ + °° duj HTiFk) + </>(u T irj) 1 ^ 

T ’ d£{k,u)/du\u= T iv k •'~°° 2?r e(k, lo + ip) ^ - F iT fc 

and 

B t) = • e^ Vkt r + °° du: (j>{^iT k ) + <Ku) ± irj) 1 c _ iujt 

T ’ de(k,uj)/duj\ UJ=Ti r k 2?r e(k,u±irj) 

III. RESULTS 

In the numerical calculations we employ an effective Skyrme potential given in Ref. [6]. 
For symmetric nuclear matter it has the form 

u( ^ A (i) +B (if- Dv2p ’ (27) 

where p(f,t) is the local nucleon density. The numerical parameters A = —356.8 MeV, 
B = +303.9 MeV, a = 1/6, and D = +130.0MeU • fm 5 are adjusted to reproduce satura¬ 
tion properties of symmetric nuclear matter: binding energy £ 0 — 15.7 MeV /nucleon, zero 
pressure at saturation density p 0 = 0.16/m -3 , compressibility modulus K = 201.0 MeV, 
and the surface energy coefficient of the Weizsacker mass formula a sur f = 18.6 MeV. Fig. 2 
shows phase diagrams corresponding to different wavelengths in temperature-density plane. 
These diagrams indicate the boundaries of spinodal unstable regions for different wave¬ 
lengths, starting from upper boundary for A = oo. The critical temperature corresponding 


(24) 

(25) 

(26) 








to this effective Skyrme potential is T c = 15.5 MeV, which occurs approximately at a density 

P = 0.3p 0 - 

We pick two reference states with initial densities p = 0.2 p 0 and p = 0.4p 0 , and calculate 
equal time correlation function of density fluctuations at two different temperatures T = 
l.OMeVand T = 5.0 MeV. In Figs.(3-6), we plot spectral intensity a(k,t) of correlation 
functions as a function of wave number k at two different times t = 30 fm/c and t = 
50fm/c for two different temperatures and initial densities as indicated above. At each 
initial density and temperature, the upper limit of the wave number k max is determined by 
the condition that the inverse growth rate of the mode vanishes, r k = 0. Dashed, dash- 
dotted and solid lines indicate the result of calculations of Eq. © with pole contributions 
only, with cut contributions only, and the total of all terms. The cut contributions include 
cut-cut contribution a cc {k,t), and the cross terms due to pole and cut parts 2a pc (k,t). From 
these figures we make two important observations: Cut terms make an important negative 
contribution during the early phase of growth, hence slowing down the growth of instabilities. 
During later times, collective poles dominate the growth of density fluctuations, and the cut 
terms representing the effects of non-collective poles do not grow in time, as discussed in 


FIG. 2. 



different wavelengths cor¬ 


responding to the potential given by Eq. (27). 
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FIG. 3. (Color online) Spectral intensity of the correlation function as a function of wave number 
at initial density p = 0.2pofm ~ 3 at time t = 30 fm/c at temperature T = 1 MeV (a) and T = 
5 MeV (b). Dotted, dashed-dotted and solid lines are results of pole, cut, and total contributions, 
respectively. 


earlier studies [1]. We also observe that both pole and cut contributions have divergent 
behavior with opposite signs, as wave numbers approach their upper limit, k —» /c max . These 
divergent behaviors cancel each other out to produce a nice regular behavior of the spectral 
intensity as a function of wave number. 

Figs. (7-8) show the equal time correlation functions a(\r — r'\,t) of density fluctuations 
as a function of distance between two space locations x = \r — r'\ at two different initial 
densities and temperatures corresponding to Figs. (3-6). The density correlation functions 
calculated with full spectral intensities are shown by solid lines in Figs. (3-6), and the results 
are plotted in Figs.(7-8) at four different times t = 20, 30,40, 50fm/c. Growth of density 
fluctuations is clearly visible in these figures. In addition, we can obtain information about 
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FIG. 4. (Color online) Spectral intensity of the correlation function as a function of wave number 
at initial density p = 0.2pofm ~ 3 at time t = 50fm/c at temperature T = 1 MeV (a) and T = 
5 MeV (b). Dotted, dashed-dotted and solid lines are results of pole, cut, and total contributions, 
respectively. 


typical size of condensation regions during the growth of fluctuations. For this purpose, we 
introduce correlation length xc as the width of the correlation function at half maximum. 
The correlation length provides a measure for the size of condensing droplets during the 
growth of fluctuations. In the correlation volume AVer = Attx^/3, the variance of local 
density fluctuations at time t is approximately given by a(xc,t). The number of nucleons 
in each correlation volume fluctuates with dispersions A Ac = A Vc^jer^xc-, t). As a result, 
the nucleon number in each correlation volume fluctuates approximately within the range 
AAo — A Ac < AA < AAo + A Ac, where AAo = A Vcp denotes the number of nucleons 
at the initial uniform state. From this analysis we notice that the spinodal decomposition 
mechanism does not necessarily lead to similar size fragmentation of the system. Rather 
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FIG. 5. (Color online) Spectral intensity of the correlation function as a function of wave number 
at initial density p = OApofm ~ 3 at time t = 30 fm/c at temperature T = 1 MeV (a) and T = 
5 MeV (b). Dotted, dashed-dotted and solid lines are results of pole, cut, and total contributions, 
respectively. 


this mechanism leads to a mixture of various cluster sizes . As seen from Figs. (7-8), 
correlation length is not very sensitive to the temperature and time evolution, but depends 
on the initial density. The correlation lengths are xc = 2.0 fm and xc = 3.0 fm for the 
initial densities p = 0.2 p 0 fm~ 3 and p = 0Ap o fm ~ 3 , respectively. As an example, we note 
that in the lower panel of Fig|8j at the initial density p = 0.4p 0 /^~ 3 and temperature 
T = 5.0Meld,the magnitude of dispersion of density fluctuations at time t = 50 fm/c is 
about yja(xc, t) = 0.04 fm~ 3 . Consequently, the number of nucleons in the correlation 
volume approximately fluctuates in the range 3 < A A < 12. 

We like to clarify that the linear response treatment provides a description for the early 
evolution of the condensation mechanism inside the spinodal region. Since the density fluc- 
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FIG. 6. (Color online) Spectral intensity of the correlation function as a function of wave number 
at initial density p = 0Apofm~ 3 at time t = 50 fm/c at temperature T = lMeV(a) and T = 
5MeV(b). Dotted, dashed-dotted and solid lines are results of pole, cut, and total contributions, 
respectively. 


tuations keep growing nearly exponentially in time, the linear response treatment of fluctua¬ 
tions beyond a limiting time does not provide a realistic description of instabilities. We can 
constraint the linear response regime approximately by the condition that the dispersion on 
nucleon density fluctuations remains below the average nucleon density in the correlation 
volume, i.e., \Jcr{xc,t) < p. As seen in Fig. (8), this condition is well satisfied for the initial 
density p = 0Ap o fm~ 3 for times up to t — 50 fm/c. In a similar manner, we can estimate 
nucleon number fluctuations in the correlation volume for the initial density p = 0.2pofm~ 3 . 
However, we note that according to the constraint, as seen in Fig. (7), at this low initial 
density the linear response regime is already exceeded beyond time t = 25 fm/c. We also 
note that, in the linear response treatment, even though density fluctuations and hence early 
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FIG. 7. (Color online) Density correlation function as a function of distance between two space 
location x = \r — r 1 | for initial density p = 0.2pofm~ 5 and temperature T = lMeV(a) and 
T = 5MeV(b) at times t = 20,30,40, 50/m/c. 


size of condensing fragments grow in the correlation volume, the matter is not expanding, 
and hence size of the correlation volume remains nearly constant in time. In the linear re¬ 
sponse description, we can approximately specify the range of early cluster size distribution, 
however we cannot identify a favored cluster size. For more realistic treatment of fragmen¬ 
tation mechanism, we need to study nonlinear evolution of the density fluctuations in the 
framework of the SMF approach. In such a non-linear evolution, it is expected that the 
exponential growth of correlation function of density fluctuations saturates. Furthermore, 
since in such a non-linear evolution the matter is expanding, the size of correlation volume 
become enlarged in time and the fragment size distribution is modified by recombination 
and breakup of clusters produced during the early phase of the evolution. But this study is 
beyond the description of the present work. 
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FIG. 8. (Color online) Density correlation function as a function of distance between two space 
location x = \f — r '\ for initial density p = 0.4po/m _3 and temperature T = lMeF(a) and T = 
5MeV(b) at times t = 20,30,40,50 fm/c. 


In previous studies [17,19-22], in the calculations of the density correlation functions we 


included only the contributions due to the collective pole in the residue integral in Eq. (10). 
As seen from dashed lines in Figs. (3-6), we encounter the problem that both pole and cut 
contributions have divergent behavior as wave numbers approach its upper limit, k —>• k max . 
To cure this problem in those calculations, we introduced a cut-off in the integral over the 


wave numbers in Eq. (18). At initial density of p = 0.2po/m 3 the cut-off value was taken 
at the local minimum point of the spectral intensity. 


In Fig. ([9]) calculations with cut-off are compared with the complete calculations of the 
density correlation function at temperatures T = lMeVand T = 5 MeV. Indeed, we see 
that the cut-off provides a good approximation at this initial density. At the initial density 
of p — OAp 0 fm ~ 3 , there is no visible local minimum of the spectral intensity. I 11 this case, we 
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FIG. 9. (Color online) Density correlation function as a function of distance between two space 
location x = \r — r '\ for initial density/? = 0.2pofm~ 3 and temperature T = 1 MeV (a) and 
T = 5MeV(b) at times t = 30,50 fm/c. Solid and dotted lines are complete calculations and pole 
contributions, respectively. 

cut the integrations over the wave number at k cut = 0.8 fmT 1 . In Fig.(10), we compare the 
density correlation function calculated with the complete spectral intensity and with those 
calculated with pole contribution with cut-off at temperatures T = lMeDand T = hMeV. 
We observe that the results obtained with pole approximations are rather sensitive to the 
cut-off wave number. 

IV. CONCLUSIONS 

We investigate early growth of density fluctuations in nuclear matter in using the SMF 
approach. In the linear response framework, employing the method of one-sided Fourier 
transform, it is possible to carry out a nearly analytical treatment of the correlation function 
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FIG. 10. (Color online) Density correlation function as a function of distance between two space 
location x = \f — r '\ for initial density p = 0.4po/m _3 and temperature T = lMeF(a) and T = 
5MeV(b) at times t = 30,50 fm/c. Solid and dotted lines are complete calculations and pole 
contributions, respectively. 


of density fluctuations. The density correlation function provides a very useful information 
about the early condensation mechanism during the liquid-gas phase transformation of the 
system. The one-sided Fourier transform method involves a contour integration in the 
complex frequency plane. In earlier investigations, this contour integral was evaluated by 
keeping only effects of the unstable collective poles, and the effects of non-collective poles 
were ignored. In the spectral intensity of the density correlation function, the collective 
pole contribution has an important drawback that it leads to a divergent behavior as the 
wave numbers approaches its upper limit, k —> /c max . In previous studies to cure this 
problem, the integration over wave numbers was terminated at a cut-off wave number. 
Since short wavelength fluctuations cannot grow in time, cut-off provides a relatively good 
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approximation. However, the complete description is always desirable. In this work, we 
consider symmetric nuclear matter in non-rclativistic framework of the SMF approach. We 
calculate the spectral intensity of the correlation function, including the effects of collective 
poles and non-collective poles in terms of the cut contribution in the complex frequency 
plane. The cut contribution also has a divergent behavior with opposite sign, as wave 
numbers approach the upper limit, k —> /c max . As a result divergent behaviors of pole and 
cut contributions cancel each other out to produce a nice regular behavior of the spectral 
intensity as a function of wave number. This allows us to have a complete description 
of the correlation function of density fluctuations. It will be very interesting to carry out 
similar calculations of correlation functions for charge asymmetric nuclear matter in the semi- 
classical and the quantal framework by employing non-relativistic and relativistic versions 
of the SMF approach. 


ACKNOWLEDGMENTS 

S.A. gratefully acknowledges TUBITAK and the Middle East Technical University for 
partial support and warm hospitality extended to him during his visits. This work is sup¬ 
ported in part by US DOE Grant No. DE-FG05-89ER40530, and in part by TUBITAK 
Grant No. 114F151. 


[1] Ph. Chornaz, M. Colonna, and J. Randrup, Phys. Rep. 389 (2004) 263. 

[2] L. Beaulieau et al., Phys. Rev. Lett. 84 (2000) 5971. 

[3] B. Borderi et al., Phys. Rev. Lett. 86 (2001) 3252. 

[4] S. Ayik and C. Gregoire, Phys. Lett. B212 (1988) 269; Nucl.Phys. A513 (1990) 187. 

[5] J. Randrupand B. Remaud, Nucl. Phys. A514 (1990)339. 

[6] V. Baran, M. Colonna, M. Di Tora, and A. B. Larionov, Nucl. Phys. A632 (1998) 287. 

[7] Y. Abe, S. Ayik, P.-G.Reinhard, and E. Suraud, Phys. Rep. 275 (1996)49. 

[8] M. Tohyama, and A. S. Umar, Phys. Lett. B549 (2002)72. 

[9] D. Lacroix, S. Ayik, and Ph. Chornaz, Prog. Part. Nucl.Phys. 52 (2004)497. 

[10] S. Ayik, Phys. Lett. B658 (2008)174. 


18 



[11] D. Lacroix, S. Ayik, and B. Yilmaz, Phys. Rev. C85 (2012) 041602(R). 

[12] D. Lacroix, D. Gambacurta, and S. Ayik, Phys. Rev. C87 (2013) 061302(R). 

[13] K.Washiyama, S. Ayik, and D. Lacroix, Phys. Rev. C80 (2009)031602(R). 

[14] B. Yilmaz, S. Ayik and D. Lacroix, Phys. Rev. C83 (2011) 064561. 

[15] S. Ayik, K. Washiyama, and D. Lacroix, Phys. Rev. C79 (2009)054606. 

[16] D. Lacroix and S. Ayik, Eur. Phys. J. A50, (2014) 95. 

[17] S. Ayik, N. Er, O. Yilmaz, and A. Gokalp , Nucl. Phys.A 812 (2008) 44. 

[18] D. Vretenar, A. V. Afanasjev, G. A. Lalazissis, and P. Ring, Phys. Rep. 409 (2005) 101. 

[19] S. Ayik, O, Yilmaz, N. Er, A. Gokalp, and P. Ring, Phys. Rev. C 80 (2009)034613. 

[20] S. Ayik, O. Yilmaz, F. Acar, B. Danisman, N. Er, and A. Gokalp, Nucl. Phys.A 859 (2011) 
73. 

[21] O. Yilmaz, S. Ayik, and A. Gokalp, Eur. Phys. J.A 47 (2011) 123. 

[22] O. Yilmaz, S. Ayik, F. Acar, S. Saatci, and A. Gokalp, Eur. Phys. J. A 49(2013) 33. 

[23] P. Bozek, Phys. Let. B383 (1996) 121. 

[24] E. M. Lifshitz and L. P. Pitaevskii, Physical Kinetics, Pergamon (1981). 

[25] H. Heiselberg, C. J. Pethickand, and D. G. Ravenhall, Phys. Rev. Lett. 61 (1988) 818. 

[26] S. Ayik, Ph. Chomaz, M. Colonna, and J. Randrup, Z. Phys.A 335 (1996) 407. 


19 



